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Abstract. The Mauritius Radio Telescope (MRT) images show systemat- 
ics in the positional errors of sources when compared to source positions in the 
Molonglo Reference Catalogue (MRC). We have applied two-dimensional homog- 
raphy to correct positional errors in the image domain and avoid re-processing 
the visibility data. Positions of bright (above 15-cr) sources, common to MRT 
and MRC catalogues, are used to set up an over-determined system to solve 
for the 2-D homography matrix. After correction, the errors are found to be 
within 10% of the beamwidth for these bright sources and the systematics are 
eliminated from the images. 

1. Introduction 

Images of about a steradian (18^ < a < 24^^30™, - 70° < 6 < - 10°) of the 
southern sky, at 151.5 MHz, using MRTludava Shankar et al.ll2002l ). were pro- 



duced bv lPandev and Udava Shankail (|2005l ) . The deconvolved images, with an 
angular resolution of 4' X 4'. 6 sec (5-|-20°.14) and rms noise of about 300 rn Jy beam~^, 
and a source catalogue of 2782 sources were published bv lPandevI (l2006l ). 



Systematics in positiona ,l errors were foun d when the positions of sources 
common to MRT and MRC ( Large et al.lll98ll ) catalogues were compared. The 



positional errors in right ascension (RA) show no systematics as a function of RA 
and the errors are within 15% of the RA beamwidth. The positional errors in 
declination (DEC) also show no systematics as a function of RA. However, errors 
in DEC show a linear gradient as a function of DEC (Figure [Ka)) and reach 
about 50% of the MRT beamwidtli in declination. Histogram in Figure [IJb) 
shows the distribution of err ors in DEC. Du e to space limitations similar plots 
for RA errors are not shown. Pandev ( 20061 ) used one-dimensional (1-D) robust 



least squares fit, estimated the errors and corrected them in the source catalog. 
However, the errors remained in the images, which impede usefulness of MRT 
images for multi-wavelength analysis of sources. 

This paper describes the application of two-dimensional (2-D) homography, 
a technique ubiquitous in the computer vision and graphics community, to cor- 
rect the errors in the image domain. At MRT, the visibility data is processed 
through several complex stages of data reduction specific to tii e array, especially, 
arising due to its non-coplanarity ( Udava Shankar et al.|[2b02l ). It was therefore 



decided to correct for the errors in the image domain and avoid re-processing 
the visibility data. Homography is used to estimate a transformation matrix 
(which includes rotation, translation and non-isotropic scaling) that accounts 
for positional errors in the linearly gridded 2-D images. 
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Figure 1. Positional error analysis, (a) Plot in the first row shows DEC 
errors against RA; no systematics are observed. Second row plot shows DEC 
errors against DEC; errors show a linear gradient as a function of DEC. (b) 
Histogram of positional errors in DEC. 



2. Homography-based correction scheme 



2-D planar homography is a non-singular linear relationship between points on 
planes. Given two sets of corresponding image points in proi ective coordinates, 
p, G P^, homography maps each p to the corresponding p' (,Hartley and Zisserman 
I2OOCII ). The homography sought is a non-singular 3x3 matrix H such that: 
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Where, and {x',y') represent (RA, DEC) of MRT and MRC sources, re- 

spectively. H is defined up to a scale factor, hence the problem has 8 degrees 
of freedom. Each point provides two equations, so a minimum of 4 point cor- 
respondences are necessary to fully constrain H. The equation contributed by 
each corresponding pair can be re-arranged to obtain: 

X (/131X + /i32y + /las) = hux + hi2y + /113 

y'ihsiX + hz2y + /iss) = h2lX + /l22y + /l23- 

A set of n such equation pairs, contributed by n point correspondences, form an 
over-determined linear system Ah = b, as shown in Eq. [3j 
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Figure 2. Block schematic of the correction scheme. Rectangular boxes 
represent processes, rounded boxes data or results. 
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Figure 3. Positional error analysis after homography-based correction, (a) 
First and second rows show DEC errors against RA and DEC, respectively, 
(b) Histogram of DEC positional errors. 



This system can be solved by least squares based estimators. To obtain 
a good estimate of the transform ation matrix we adopted the normalization 
scheme proposed by [Hartley ( 19971 ). This ensures freedom on arbitrary choices of 
scale and coordinate origin, leading to algebraic minimization in a fixed canonical 
frame. 

Figure [2] shows the block schematic of the correction scheme. We identify 
bright (above 15-a) sources common to the MRT and MRC catalogues. There 
are about 400 common sources in the steradian of the sky under consideration. 
A single 3x3 homography matrix (refer equation [T|) was estimated from this 
population. Each pixel from the images was then projected to a new position 
using this matrix. 



3. Results 



Figure [3] shows positional errors in DEC after homography-based correction. 
Compare these plots with Figure [TJ Clearly, homography has managed to remove 
systematics and reduce the errors to well within 10% of the beamwidth for 
sources above 15-a. 

Figures (D^a) and Hl^b) show MRT contours before and after correction, re- 
spectively, overlaid on SUMSS (Sydney University Molonglo Sky Survey) image 
(|Mauch et al.ll2003l 1. for a source around -67° DEC. The corrected MRT image 
contours in Figure |3|^b) overlap well with the source in SUMSS image. 
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Figure 4. MRT contours overlaid on SUMSS image, (a) MRT contours 
before correction and, (b) MRT contours after homography-based correction. 

4. Conclusions 

The homography-based correction was able to correct for systematics in posi- 
tional errors in the image domain and the errors are now well within 10% of the 
beamwidth for sources above 15-0". The corrected MRT images are available for 
download at http://www.rri.res.in/surveys/MRT, 

Positional error analysis clearly shows that MRT images are stretched in 
declination (about 1 part in 1000). This translates to a compression of the 
baseline scale, in the visibility domain, which cued towards possible errors in 
our estimation of the array geometry. By formulating a simple linear system, 
using instrumental phases estimated from three well separated c alibrators whose 



posit ions are well known, the array geometry was re-estimated (jThomson et al 
I2OOII ). The estimates show an error of about 1 mm/m. In other words, the error 
is about half a wavelength at 150 MHz (1 m) for a 1 km baseline. This matches 
with the observed stretching of MRT images. 
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